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The exclusive diffractive production of vector mesons and real photons in ep col- 
lisions has been studied at HERA in a wide kinematic range. Here we present the 
most recent experimental results together with a Regge-type model. We deduce the 
Pomeranchuk trajectory (Pomeron) by analyzing the HERA data on deeply virtual 
D ■ Compton scattering (DVCS), and then discuss its basic properties, namely its appar- 

O ■ ent "hardness" and its "non-flat" behavior, different from the claims of some authors. 
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1 Introduction 

The diffractive scattering is a process where the colliding particles scatter at very small 
angles and without any color flux in the final state. This involves a propagator carrying 
the vacuum quantum numbers, called Pomeron and described, in the soft regime, within the 
Regge theory. Since the first operation period in 1992, ZEUS and HI, the two experiments 
dedicated to the DIS physics at HERA, observed that a big amount (~ 10%) of lepton-proton 
DIS events had a diffractive origin opening a new area of studies in diffractive production 
mechanism, providing a hard scale which can be varied over a wide range and therefore it is 
an ideal testing for QCD models. 

In particular, the diffractive production of Vector Mesons (VMs) and real photons in ep 
collisions allows to study the transition from the soft to the hard regime in strong inter- 
actions. The hard regime (high energy and low Bjorken-x) is dominated by the exchange 
of a hard Pomeron sensitive to the gluon content and well described by perturbative QCD 
(pQCD), while in the soft regime (low- re) the interaction is well described within the Regge 
phenomenology. Indicating with Q 2 the virtuality of the exchanged photon and with M 2 the 
square mass of the produced VM, HERA data suggested a universal hard scale, Q 2 + M 2 , for 
the diffractive exclusive pruduction of VMs and real photons, which indicates the transition 
from the soft to the hard regime. 



2 Q 2 and W dependence of the cross section 

A new precision measurement of the reaction 7*p — > p°p was published by ZEUS [lj. It 
was found that the cross section falls steeply with the increasing of Q 2 but, unlike it was 
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observed for the J / ip electroproduction [2], [3] , it cannot be described by a simple propagator 
term like a oc (Q 2 + M 2 )~ n , in particular an n value increasing with Q 2 appears to be 
favored. Figure [T] reports the cross section for the p° electroproduction versus Q 2 compared 
with several theoretical predictions: the KWM model [I] based on the saturation model, the 
FSS model [5] with and without saturation and the DF model [6]. None of the available 
models gives a good description of the data over the full kinematic range of the measurement. 



ZEUS 




Q 2 (GeV 2 ) 

Figure 1: The j*p — ► p°p cross section as a function of Q 2 measured at W = 90 GeV 2 and 
comared in (a) and (b) with different models as described in the text. 

The soft to hard transition can be observed looking at the dependence of the VMs 
photoproduction (Q 2 = 0) cross section from the j*p centre of mass energy, W, where the 
scale is provided by M 2 . Figure [2] collects the <r(7*p — ► Vp) as a function of W from the 
lightest vector meson, p°, to the heaviest, T, compared to the total cross section. The cross 
section rises with the energy as W , where the 6 exponent increases with the hard scale 
M 2 as expected for a transition from the soft to the hard regime. New results on the T 
photoproduction [7], recently published by ZEUS, confirmed the steeper rise of a(W) for 
higher vector meson masses. 

The transition from the soft to the hard regime can also be studied varying Q 2 . Recent 
results were achieved by HI [?] and ZEUS |9j for the exclusive production of a real photon, 
the Deeply Virtual Compton Scattering (DVCS), where the hard scale is provided only by 
the photon virtuality, Q 2 . Figure [3] shows the HI (left) and the ZEUS (right) results. The 
steep rise with W of the cross section even at low-Q 2 , seems to suggest that the most sensitive 
part to the soft scale comes from the wave function of the pruduced VM. A similar result 
was obtained for the J/ip electroproduction [21 E]- 

The electroproduction of a large variety of VMs was studied at different Q 2 values and 
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Figure 2: The W dependence of the cross section for exclusive VM photoproduction together 
with the total photoproduction cross section. Lines are the result of a W s fit to the data at 
high PU-energy values. 
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Figure 3: The W dependence of the cross section for a DVCS process. Lines come from a W s 
fit to the data. Left: the HI measurement of the 5 slope as a function of Q 2 . Right: the new 
ZEUS preliminary measurement at low Q 2 (dots) together with the published measurements 
(squares). 
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the corresponding slope 5 is reported in Fig. H] (left) versus the scale Q 2 + M 2 , including the 
DVCS measurements. Data show a logarithmic shape 5 oc ln(Q 2 + M 2 ) and the behaviour 
seems to be universal with 5 increasing from 0.2 at low scale, as expected from a soft Pomeron 
exchange [10] to ~ 0.8 at large scale values. 
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Figure 4: The dependence on the hard scale Q 2 + M 2 of the value 5 (left) extracted from a 
fit W s and of the slope B (b in the figure lable) (right) extracted from a fit ^ oc e B '*' for the 
exclusive VM electroproduction. DVCS is also included. 



3 t dependence of the cross section 

The differential cross section as a function of t, the four-momentum tranfer at the proton 
vertex, can be parametrised by an exponential fit: 4|r cx e s '*L Figure H] (right) reports the 
collection of the B values versus the scale Q 2 + M 2 for the electroproduction of VMs and 
DVCS, with B decreasing from ~ 11 GeV~ 2 to ~ 5 GeV~ 2 as expected in hard regime. 

The measurement of da/d\t\ for the DVCS process, recentrly published by the HI Col- 
lab [?], where t was obtained from the transverse momentum distribution of the photon, 
studied B versus Q 2 and W as shown in Fig. [51 B seems to decrease with Q 2 up to 
the value expected for a hard process but it doesn't depend on W. A new preliminary 
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Figure 5: The t slope parameter B {b in the figure lable) as a function of Q 2 (left) and W 
(right). 
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ZEUS measurement [9] of da/d\t\ has been achieved from a direct measurement of the pro- 
ton final state of using a spectrometer based on the roman pot thechnique. The result 
B = 4.4 ± 1.3 (stat.) ± 0.4 (syst.) GeV~ 2 } measured at Q 2 = 5.2 GeV 2 and W = 104 GeV, 
is consistent, within the large uncertainties due to the low acceptance of the spectrometer, 
with the HI result [?] of B = 5.45 ± 0.19 (stat.) ± 0.34 (syst.) GeV" 2 at Q 2 = 8 GeV 2 and 
W = 82 GeV. 

Since the B value can be related via a Fourier transform to the impact parameter and 
assuming that the exclusive process in the hard regime is dominated by gluons, the relation 
(r 2 ) = 2b(hc) 2 can be used to obtain the radius of the gluon confinement area in the proton. 
b ~ 5 GeV 2 corresponds to (r 2 ) ~ 0.6 fm smaller than the proton radius (~ 0.8 fm) 
indicating that the gluons are well contained within the charge-radius of the proton. 



Exclusive production of vector states (real photon, VMs, lepton pairs) via deeply virtual 
scattering, ep — > epV, is interesting for many reasons, above all as a source in extracting 
informations about General Parton Distributions (GPDs). Of particular interest is DVCS in 
which the outgoing real photon, interferring with the photon coming from bremsstrahlung 
(Bethe-Heitler process), offers a holografic picture of the nucleon. 

A Q 2 -dependence of generalized parton distributions (GPD) can be found from QCD 
evolution, similar, although less explored than that of ordinary parton distribution, obeying 
the DGLAP evolution equation Q. 

In the present paper, we analize and test an explicite model for DVCS. The model can 
be used to infer the GPD by a relevant deconvolution procedure. 

After the shut-down of the electron-proton collider HERA, the HI and ZEUS collabora- 
tions have left a huge heritage of data still waiting for a better understanding. In particu- 
lar, it concerns the properties of the Pomeron trajectory as seen in VM electroproduction, 
ep — > eVp, and in DVCS, ep — > e^p. There are many papers discussing in details the form 
and the values of the parameters of the Pomeron trajectory as well as their possible Q 2 depen- 
dence (for recent reviews see, e.g. [llj). In general they introduce a linear Regge-trajectory 
of the Pomeron 



where t is the squared four-momentum transferred at the proton vertex and Q 2 is the vir- 
tually of the exchanged photon. For the "effective" Pomeron they uses the standard Regge 
pole parametrization of the scattering amplitude 



where s = W 2 is the squared Yp centre-of-mass energy and B(t, Q 2 ) is related to the radius 
associated with the proton vertex. 

The problem of constructing Regge-type models for currents (lepton-hadron scattering) 
or the off-mass- shall continuation of the analytic ^-matrix has already a long history, still 
lacking a consistent solution. The concepts of the analytic S"-matrix formally may even be 
incompatible with, or inapplicable to processes like Compton scattering. Nevertheless, the 
need for a theoretical framework to describe high-energy lepton-hadron scattering served 

1 for a discussion of GPD evolution see the Appendix 



4 The Pomeron in DVCS at HERA 



a P (t,Q 2 ) = a (Q 2 ) + a'(Q 2 )t, 



(1) 
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as an "excuse" in using Q 2 — dependent Regge-type models for small x (large s ~ Q 2 /x) 
processes, like Compton scattering. Moreover, to meet the apparent acceleration of the rise of 
the structure functions with 1/x towards larger Q 2 , a Q 2 dependent Pomeron was introduced, 
in complete discords with Regge factorization, in which the dependence on the mass or 
virtuality of the external particle (photon or meson in our case) can enter at most through 
the relevant vertex function. This inconsistency was circumvented by calling this object 
"effective Pomeron" (or, more generally, a Reggeon), implying that it accommodates by more 
than a single Regge exchange, without specifying the content. Very often, however, this was 
ignored and the parameters of the Pomeron have been extracted from simple formula like in 
Eq. [TJ Most of the conclusions, especially the appearance a large and ever increasing intercept 
(hard Pomeron) and its small slope (flat Pomeron) were used to confirm perturbative QCD. 

At the same time, it was shown [12] that the inclusion of a sub-leading Regge contribution 
modifies the parameters of the Pomeron, in particular the observed rise of the structure 
function partly is due to the decrease of the sub-leading contribution, and the Q 2 dependent 
residue functions can provide for the apparent "hardening" of the Pomeron. 

Recently a number of Regge /Pomeron-type models for DVCS appeared in the litera- 
ture [IHl HH1 HEl HTj. 

In the present paper, we analyze the high-energy data on DVCS collected by the HI and 
ZEUS detectors at HERA [HI El UHl UH] - To do so, we use a model for the Pomeron elaborated 
and fitted to the data in Ref. [13] . 

In most of the Regge/pole models used in analyzing the experimental data, a linear 
Pomeron trajectory is used. Although linear trajectories contradict to the postural of the 
^-matrix (analyticity and unitarity), to perturbative QCD calculations (BFKL) [14], as well 
as disagree with precize fits to the data, they can serve as a simple approximation to the 
observed phenomena at small \t\ [13] however the linear behavior is replaced by a logarithmic 
one0. 

For this reason we use a simple Regge pole model for the Pomeron, however, contrary to 
most of the known models, our Pomeron trajectory is essentially nonlinear. 
The scattering amplitude has the form 

A(s,t,Q 2 ) rp ^ lp = -A e ba ^e b ^\-zs/s r^ = - Ao e^ L ^ +b ^ z \ (3) 
where L = ln(— is/so). The trajectory at the pIPp vertex is 

a(t) = «o — Oil ln(l — a2t). (4) 

whereas the trajectory at the 7*/P7 vertex is 

(3(z) = ao — «i ln(l — ac2z), (5) 

with z = t — Q 2 - a new variable introduced in Ref. [13]. Notice that the presence of the 
"minus" sign in Eq. [3] is important for the linear forms, e.g. ImA(t = 0), proportianal to the 
the total cross section, or for the ratio p = ReA/ ImA, but it is irrelevant for the squared 
modulus we are interested in. 

Similar to [15] . we consider only the helicity conserving amplitude. For not too large Q 2 
the contribution from longitudinal photons is small (it vanishes for Q 2 = 0). Moreover, at 
the high energies typical of the HERA collider, the amplitude is dominated by the helicity 
conserving Pomeron exchange and, since the final photon is real and transverse, the initial 

2 For further arguments see the Ref. [13] 
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one is also transverse and the helicity is conserved. The electroproduction of vector mesons 
requires the account for both the longitudinal and transverse cross sections. 

Generally speaking, any DVCS reaction admits the contribution from a subleading Regge 
exchanges, most important of which is the / trajectory, made of a quark and an antiquark, 
not gluons like the Pomeron. It differs only by the values of the Regge-trajectory parameters 
(lower intercept and steeper slope). 

The cross section can be calculated as 

^M,Q 2 ) = y 2 \A(s,t,Q 2 )\ 2 . (6) 
with the slope of the da/d\t\ cross section coming from 

The model contains quite a number of parameters but some of them (sq, cti, 0.2) can be 
fixed by theoretical constraints, as explained in Ref. [13j. Fits to the DVCS data collected 
at HERA [TBI [8j [19] were performed and the results are shown in Fig. [BJ with the values of 
the fit parameters, A , b, «o quoted in Table [HI Note that a(W) is sensitive to the Pomeron 
intercept, «o> an d the corresponding fit was performed by keeping bp = b = 1.0 fixed, while 
the cr(Q 2 ) is sensitive to b and the corresponding fit was done at fixed cto = 1.2. Figures [6K,b 
show that the model can reproduce all the features of data versus Q 2 and W = \fs. Results 
suggest that the Pomeron in DVCS has an intercept, a® ~ 1.2 graeter than that in hadronic 
reactions, but a slope, a' ~ 0.25, typical of hadronic scattering and not so flat as it was 
observed for the diffractive electroproduction of vector mesons {a' ~ 0.1 ~ O.bathh') in hard 
regime. 



Parameter 


@ DVCS vs Q 2 


&DVCS VS W 


da/d\t\ HI 


da/d\t\ ZEUS 


A 


0.13 ± 0.01 


0.10 ± 0.01 


0.17 ± 0.01 


0.44 ± 0.10 


a 


fixed at 1.2 


1.23 ± 0.006 


fixed at 1.2 


fixed at 1.2 


b 


0.96 ± 0.07 


fixed at 1.0 


fixed at 1.0 


fixed at 1.0 


x 2 


1.39 


1.21 


16.8 


0.3 



Table 1: Values of fitted parameters and the corrispondent \ 2 value. 

The fit of da/d\t\ was performed with all the parameters fixed excepted for the normalisa- 
tion. The model does not agree with the HI measurements (see Fig. [6b) but it is compatible 
with the new ZEUS preliminary results (see Fig. [6H). The t variable was calculated in HI 
by the transverse four-momenta of the scattered electron and the real photon, using the ap- 
proximation t ~ \Pt £ + Pt~,\ 2 , while in ZEUS a particular silicon microstrips spectrometer, 
based on the roman pots technique, was used in order to have a direct measurement of the 
scattered proton momentum p' and then t being calculated from the quantity: Xl = Pz , 

^ ^ " Pbeam 

by the phormula: t = Px ^ y + P 2 eam • ^ ne rea Uy l° w acceptance of this spectrometer 

is the reason of the poor statistics of the ZEUS data, howewer it offers a really pure selection 
of diffractive events, not affected by any not-diffractive background. The agreement with 
the ZEUS preliminary measurements encourages us to be still confident in our predictions, 
till they will be checked by the new HERA data analyses (now in progress with and without 
roman pots spectrometer). 



7 



c 



£10 

T 

a 



10 



1 1 1 1 1 1 1 1 1 
■ • 

1= ■ 


i 1 i i i 1 i i i 
ZEUS f \ 
HI 96/00 


- ^ A 


HI 05/06 : 




This model _ 


! Q 2 = 8GeV 2 

" 


« • \ 

■ 



20 40 



60 



80 100 

Q 2 (GeV 2 ) 



> 
O 

> 10 
e 



T 

o, i 

b 

6 
■o 



10 



1 1 1 


1 1 1 1 1 
■ 


i i 1 i i i 1 i 
HI 96/00 


_ ■ 




This model — 










8 GeV 2 




vv = 


82 GeV 




■(c) 
, , , 





! 





c 



12 



£ 10 

T 

p, 8 



> 

a 
is 

c 

? 

T 



e 
■a 



6 
4 
2 


10 : 
10 



'I 1 1 1 I 1 1 1 

• ZEUS 
■ HI 96/00 



HI 05/06 I 



This model 



(b) 




0.2 



0.4 



0.6 0.8 



10 



_L 



W = 82 GeV 





40 60 80 100 120 140 

W (GeV) 



: ' ' ' 1 

: ^ 




i 1 








» 












- Q 2 = 5.2 GeV 2 






r W = 104 GeV 






• ZEUS 00 prel. 




This model 




(d) : 



0.2 



0.4 



0.6 



Itl (GeV z 



Itl (GeV") 



Figure 6: Upper panel: Energy (a) and Q 2 (b) dependence of the DVCS cross section, the model 
is fitted to the HERA data. Lower panel: the prediction for the t dependence of the cross section is 
compared with HI (c) and ZEUS (d) data with only the normalisation kept as a free fit parameter. 
Error bars include both statistical and systematic uncertainties summed in quadrature. The error 
bars include both the statistical and systematic uncertainties summed in quadrature. 



The slope of da/d\t\, calculated according to Eq. [3, is depicted in Fig. [7] as a function of 
Q 2 and W and is compared with de HERA measurements. The local slope is predicted to 
slightly rise with W but to be almost independent on Q 2 . 

Eq. [7] was then used to make a collection of figures showing the slope dependence for 
different values of W and t. The t dependence of the local slope calculated from Eq. [7J is 
shown in Fig. [8^, at fixed Q 2 = 4 GeV 2 for three different values of energy. 

Figure [8b depicts the Q 2 dependence of the local slope at W = 60 GeV for several values 
of t. 

Figure [Hfc shows the enerrgy dependence of the local slope at Q 2 = 4 GeV 2 for several 
values of t. 
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Figure 7: Slope of the cross section vs Q 2 and W as calculated from Eq. [?| compared with the 
HERA data. 



5 Conclusions 

• After the shut-down of the HERA Collider, the ZEUS and HI experiments have left 
a huge heritage of experimental results on the diffractive production of vector mesons 
and real photons still waiting for a better understanding. Nevertheless, new analyses 
on HERA data are still in rpogress. 

• We show that data on DVCS from HERA can be fitted with a semi-hard Pomeron, 
not being flat. 

• The Pomeron trajectory could be a non-linear function; it is nearly linear at small 
and moderate values of \t\, leveling off at large \t\. The same logarithmic trajectory 
could be used in hadronic reactions, for example in analyzing, with higher precision, 
the future diffractive exclusive data from LHC, where details, such as the two-pion 
threshold [32], neglected here, should be taken into account. 

• Here we used a simple "minimal" model, with a small number of the free parameters, 
which, in view of the limited number of data points, made possible the convergence 
of the minimization procedure. A more detailed analyses should include also the indi- 
vidual contributions from the longitudinal and transverse cross sections, though their 
ratio R, spin effects, reaction ratios and comparison with possible QCD predictions, 
especially concerning the evolution in Q 2 (see next item and Appendix). 

• A complete solution for the Q 2 evolution in DVCS has not been yet found, although 
partial results are known from the literature. In a simplified, pragmatic approach the 



9 



o 

~4.5 
P3 



0.2 0.4 



W = 60 GeV 



I , i i L_ 

0.6 0.8 



1.2 



1.4 



Itl (GeV z ) 

(b) 



t = 0.0GeV 2 

t = 0.2 GeV 2 

t = 0.4GeV 2 

t = 0.6GeV 2 

___ t = 0.8GeV ! 

t= 1.0 GeV 2 

I i . i . I .... I ... i I i i i i I ii 

5 10 15 20 25 



30 



Q 2 (GeV 2 ) 




140 

W (GeV) 



Figure 8: Slope of the cross section as calculated from Eq. OJ for several values of W, Q 2 and t. 



Q 2 and t dependence in a DVCS (actually not factorizable!), reduces the problem to the 
evolution of ordinary DIS by the DGLAP procedure. These attempts are summarized 
in Appendix. 

The model and the fitting procedure presented in this paper can be extended to high- 
energy vector meson production. We intend to do so in a forthcoming publication. 
Let us notice also that any consistent extension/interpolation of the Regge approach 
to low energies (resonance region) should involve dual models, as it was shown [33] for 
the Pomeron component in case of J/ty scattering. 
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APPENDIX. GPD and QCD evolution of DVCS 

From Eq. [3] we get the real and imaginary parts of the DVCS amplitude at Q 2 = Qq, that 
will be useful for the evolution: 

ImA DVCS (x Bj ,t,Ql) = sm [7ia(t)/2]G(t,Ql)(-^-) a{t \ (8) 



ReA DVCS (x Bj ,t,Q 2 ) = -cos[n a (t)/2]G(t,Ql)(-^y {t \ (9) 

^S XBj / 

where 



G ^ = e b{*(t)+m,Qi)]. (io) 

Skewendess is defined in terms of the Bjorken variable x B j as £ ~ 2x B j — x B j and v. v., 
^■-2^/(1 + 0- 

The aim is to evolve this amplitude to higher values of Q 2 , Q 2 > Qq ~ lGev 2 and, 
at the same time, to grasp the most important properties of the corresponding generalized 
parton distributions (GPD). 

The singlet combination of GPD, corresponding to the exchange in the t— channel of a 
C — +1 charge conjugation, even signature reggeon is [27] 

H^{x, f , t) = H«(x, f , t) - H«{-x, £, t), (11) 

with H q ( + \x, 0,0) = q(x) + q(x), where q(x) is the usual quark distribution for quark q. 
Its evolution is well known. The "valence" combination, odd signature, is even under the 
change: x — *> —x. 

The point x = £, where DGLAP (|x| > £) [21 J and ERBL (\x\ < £) [22] regimes meet, is 
directly accessible to to experiment since, at leading order in a s , the imaginary amplitude 
part of the amplitude A for DVCS is proportional to a GPD at x = £ 



ImA DVCS (x Bj ,t,Q 2 ) = -7rY,e 2 q H^ +) (Z,Z,t,Ql), (12) 



q 



where H 9 ^ + '(^, £, t, Q 2 ,) is known from Eq. [HI Unfortunately, the evolution equation for 
H q{+) (£, £, t, Ql) involves H^+\x, £, t, Qg) for ar in the interval [0, 1] (for DVCS), to be known. 

An approximate solution based on conformal invariance, valid at one loop level is possi- 
ble [23]. The off- forward matrix elements of a set of twist-two operators, whose leading- log 
evolution is diagonal due to conformal symmetry, are the Gegenbauer polynomials. In the 
forward kinematics, the Gegenbauer polynomials reduce to polynomials in x B j. In Ref. [25] 
an explicit transformation has been constructed relating the off-forward and the forward 
evolutions. 

Usually it is very difficult to apply the transformation constructed in Ref. [25] to a 
realistic problem but, in our case, the simplicity of the input at makes possible an 
explicit asymptotic solution for small £ B j and £ ~ x B j/2. It has been shown in Ref. [26] that, 
in the phenomeno logically important small-£ and \t\ region, the off-diagonal distributions 
are determined unambiguously by the small-x behavior of the diagonal parton distributions. 
By setting 

H(x,0 = H(x^,t,u, 2 ), {-l<x Bj <l), (13) 
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where the values of t and £ do not change in the evolution to a higher scale /z 2 , it was shown 
in Ref. [26j that at small £ 







j dx' 






7T JO 



d UA 

dx' \ \x'\ 



y{s) 



4s(l - s) 



(14) 
(15) 



x + £(l-2s 

The proof of Eq. [TU uses the properties of the diagonal parton distributions 

q{x' -> 1) -> 0, g*(a/) = -g s (-x'), g ns (x') = q ns {-x'), 

that corresponds to the symmetry relations 

H i {x,0 = H i (x,-£), (i = (q,w),(q,ns), 

H*'(x,Z) = H*'{-x,Z), 

H q ' ns (x,0 = H q > ns {-x,0, 

where, as usual, s stays for "singlet" quark contribution and n s means "non-singlet". At small 
£, the same anomalous dimensions 7^ control both the diagonal and off-diagonal evolution. 

Consider now the prediction of [26] for the off-diagonal distribution H q (x,l;) by making 
the physically reasonable small-x assumption that 

xq(x) = N q x~ x ", (16) 

where X q may depend also on t, X q = X q (t). Then eq. (18) can be integrated 

A„+l 



T(A + 5/2) 2 
T(A + 2) ^Jo 



ds 



4s(l-s) 
x + £(l-2s) 



(17) 



where, for singlet quarks (i.e. when X q < 1) the integral is a principal value integral and 
Eq. [T7I becomes integrable for any X q < 1. 

The distribution H Q (x,^), for small £, has the form 



with 



H q (x, = r^F^x/o 



where F q (z) is a function determined by Eq.[T71 Similarly, if xg(x) = N g x~ Xq , then H g (^, £) = 
£~ As F g (l) at small £. Viceversa, since the amplitude in Eqs. El and [H] is of Regge form and 
is proportional to for small £ and fixed Q 2 , then 

y £e 2 9 H(U) = --ImA DVCS ocC 



-ait) 



and the parton distribution that appears in Eq. O has the form 

xq(x) oc 

where xq(x) can represent ^(x) if we start from Eq. [T8j 



(19) 
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An approximate analytic evolution of the input singlet parton distribution is well known if 
the input is of the Regge-like form. Since only the diagonal parton distributions are involved 
in the calculations, one can resort to libraries that generate parton distributions but, since 
we are interested in the smali-£ region. Equation [TH then determines the Q 2 dependence 
of H q (x,C,) and, finally, Eq. [18] gives the scale dependence of ImA DV cs- The knowledge of 
H q ( + '(x,£) fixes also, at leading order in a, since 

Wf, t, Q>) = X>^' n {, t, Q 2 ) + j^-j) . (20) 

Similar results were obtained in a recent paper [27] . 
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